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ABSTRACT 

We revisit the evaporation process of gas from dwarf galaxies after it has been 
photo-ionized by the UV flux from the first stars and AGNs and heated to T ~ 10 4 or 
2 x 10 4 K respectively. Earlier estimates, based on the balance between pressure and 
gravity, indicated that dark haloes of virial velocity lower than V^ V ap — 11 — 13 kms -1 
have lost most of their gas in a dynamical time. We follow the continuous evaporation 
by a thermal wind during the period when the ionizing flux was effective. We find 
that the critical virial velocity for significant evaporation is significantly higher. For 
example, if the ionization starts at z lon = 10 and is maintained until z = 2, a mass loss 
of one e-fold occurs in haloes of Kvap — 25 (or 35 kms -1 ) for T ~ 10 4 K (or 2 x 10 4 K). 
Haloes of V CV ap — 21kms -1 (or 29 kms -1 ) lose one e-fold within the first Hubble 
time at z = 10. Any dwarf galaxies with virial velocities smaller than V cvgup must have 
formed their stars from a small fraction of their gas before Zi on , and then lost the rest 
of the gas by photo-evaporation. This may explain the gas-poor, low surface brightness 
dwarf spheroidal galaxies. By z < 1, most of the IGM gas was evaporated at least once 
form dwarf galaxies, thus providing a lower bound to its metallicity. 
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1 INTRODUCTION 

The first stars and AGNs generate a flux of UV radiation 
which photoionizes most o f the hydrogen in the u niverse by 
zion ~ 7 (see a review bv lBarkana fc Loebl 1200 ll) . The gas 
is heated to a temperature of Tj on ~ 10 4 K and is kept at 
this temperature as long as the ionizing flux persists. The 
actual temperature mostly depends on the hardness of the 
radiation flux and therefore on whether its origin is in stars 
or AGNs, respectively, but also on the flux itself. The total 
emissivity of the observed quasars at z — 3 is probably suffi- 
cient to have completely reionized the low density intergalac- 
tic Helium already before that redsh ift, thereby increasing 
the temperature to Ti o n ~ 2 x 10 4 K (iMiralda-EscuddfeOOOl; 
iBarkana fc Loebll200l() . Dark-matter haloes with an escape 
velocity lower than the associated thermal velocity of the 
ionized gas, on the order of 10 kms -1 , cannot retain this hot 
gas — it would photo-evaporate. IBarkana fc Loebl (I1999T) 
have estimated the gas loss based on the balance between 
pressure and gravity, and obtained that the haloes that lose 
50% of their gas in a dynamical time are of a typical virial 
velocity of V C va P — 11 — 13 kms -1 . 

Based on this low estimate of the critical velocity for 
photo-evaporation, one could have concluded that photo- 
evaporation is not a very important process in galactic his- 
tory. This is because a similar lower bound on the virial ve- 



locity of haloes that can form stars is obtained from simple 
cooling arguments. The radiative cooling rate, which is very 
high just above 10 4 K due to atomic recombination, becomes 
extremely low just below 10 4 K. The only available cooling 
agents in this range are H2 molecules, which are both inef- 
ficient and fragile in the pre sence of even a weak UV flux 
llHaiman. Rees fc Loebl IT996h This implies that haloes of 
virial velocity below ~ 10 kms -1 cannot have their gas cool 
to form stars independently of the radiative feedback and 
the associated reionization. Indeed, there is an indication 
for such a lower b ound in the popula tion of Local Group 
dwarf galaxies fsee lWoo fc Dekell2003f) . 

Another scale associated with the radiative feedback 
process is the Jeans scale. Once the IGM is photo-ionized, 
its pressure shuts off any further gas infall into haloes of 
virial velocit ies up to ~ 30 kms -1 (see a review of numeri- 
cal results in lBarkana fc LoebfeOOll §6.5). Gas could resume 
falling into small halos after z ~ 1 — 2 when the UV back- 
ground flux declined sufficiently llBabul fc R ccs 1992), but 
only halos of V > 20 — 25 kms -1 can form mole cular hy- 
drogen by z ~ 1 llKepner. Babul fc Sperge1ll997l) . This is 
therefore a lower bound for galaxies unless they managed to 
have their gas fall into virialized haloes before Zi on . 

The observed internal velocities in Local-Group dwarf 
galaxies indicate that dwarf galaxies with virial velocities 
in the range 10 — 30 kms -1 do exist. However, a compari- 
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son of the galaxy luminosity function at the faint end (with 
a logarithmic slope close to —1) and the halo mass func- 
tion predicted by the standard ACDM cosmological model 
(with a log slope near —2) indicates that some haloes in 
this range must have remained dark with no stars in them. 
Another relevant observation is that most of the galaxies in 
this range are dwarf spheroidals, showing only little gas and 
current star formation. This calls for a more careful analysis 
of the fate of the ionized gas in such haloes. Since an effec- 
tive i onizing flux may persist until z ~ 1 — 2 feabul fc Reesl 
1992), the evaporation process is likely to continue for many 
dynamical times in the form of a continuous thermal wind, 
and thus end up with more mass loss than estimated by 
iBarkana fc Loea 1^99), and correspondingly with a higher 
critical velocity \4 vap for significant evaporation. This wind 
is analogous to the Parker wind from the sun, which is re- 
sponsible f or a mass loss of 10~ 14 Mpj yr and is well under- 
stood fe.g-. lLamers fc Cassinellilll999l) . With a higher V C va P , 
the photo-evaporation would become a major player in the 
formation of dwarf galaxies. In haloes smaller than this crit- 
ical velocity (but larger than ~ lOkms -1 ), stars can form 
only before Zi on , and any left-over gas cannot be retained. 
This could lead to the formation of gas-poor dwarfs, not un- 
like the dwarf spheroidals observed in the Local Group, or 
to completely dark haloes also in the range 10 < V < V evs , p . 
If a significant fraction of the universal gas was able to evap- 
orate from haloes that managed to form stars earlier, then 
the evaporation process might also have had implications on 
the metal enrichment of the IGM. We therefore attempt to 
evaluate the critical virial velocity for evaporation by pur- 
suing a dynamical calculation of the evaporation process. 

In f]5]we present the basic treatment of thermal winds 
applicable to non-point like mass distributions. In SjUJwe ap- 
ply the wind analysis to cosmological galactic haloes. In [2] 
we discuss our results. 



2 WIND ANALYSIS 

We wish to compute the mass-loss rate in a steady 
isothermal thermal wind in a galactic halo. The formu- 
lation is inspired by the sta ndard problem of isother- 
mal winds from stellar ob jects iParkedll96(J : lBrandd|l97Ct 
lLamers ~ Cassinclli Il999t) . The difference is largely in the 
nomenclature. The main mathematical modification that we 
have to consider is that the haloes are not point like ob- 
jects, that is, the gravitational force per unit mass is now 
—GM(r)/r 2 i with a monotonically increasing M(r). 

We begin by considering an arbitrary spherical potential 
4>(r). The gas is assumed to be heat ed and kept at a constant 
temperature T, typically ~ 10 4 K. iBarkana fc Loebl (ll99St) 
calculated the amount of gas that unbinds "instantaneously" 
as a result of the extra thermal energy. All the leftover gas, 
which energetically cannot escape, was assumed to settle 
dow n in the haloe s' potential well. 

iParkerl (Il96(t> found that even if the gas is energeti- 
cally bound, a wind can develop. One way of looking at 
it is that irrespective of how deep the potential well is, an 
exponentially small fraction of the gas can escape. Without 
continuous heating, a thermal wind would progressively cool 
as it expands, the pressure pushing it out would decrease, 
and the wind would die out. On the other hand, if the gas 



can constantly be heated to compensate for the adiabatic 
cooling, it is easier for the gas to be expelled, and a large 
flux can be obtained and maintained. 

We assume in the following analysis that the wind is 
kept at a fixed temperature. This follows from the assump- 
tion that over the time scale it takes the gas to accelerate 
out of the potential well, the radiative equilibrium with the 
ionizing radiation can keep the gas warm. We will check the 
consistency of our results with this assumption and discuss 
its meaning in H3.4I 

The basic equations are the continuity equation, the 
Euler equation of motion, and the equation of state of an 
isothermal gas: 



dp 
dt 



= -V-(pv), 



flv 

p-jjj- = -VP+ fgrav, 
P = CsP, 



(1) 

(2) 
(3) 



where c s is the isothermal speed of sound. Under the as- 
sumption of spherical symmetry, the equations become 



4(rV)'. 



/ 2 P 

V + VV = —C s — 



(4) 
(5) 



where <j>(r) is the gravitational potential normalized to at 
r = 0: 



(r) 



GM{r) 



dr. 



(6) 



and where dot and prime denote Eulerian partial derivatives 
with respect to t and r respectively. 

We now assume steady state, where M is the same at 
every r. Writing M = J 4nr 2 pdr, and using continuity, 
we obtain M = — A-kt 2 pv. Steady state therefore implies 
r 2 pv =const., and then from continuity p = 0. It also re- 
quires v = 0. The equations then combine to the wind equa- 
tion: 



2ci 



^-)v' = -<j>'(r) + 



(7) 



Any steady state theory with a transonic solution, such 
as a thermal wind, necessarily requires that the transition 
between the sub-sonic and super-sonic flow will either take 
place at a shock (which is of no interest for us here) or that 
the sonic point will coincide with the critical point where the 
net radial forces vanish. This is because a transonic solution 
to eq. |7| implies that the left-hand side changes its sign, 
and therefore the right-hand side has to change sign at the 
same radius. 1 Thus, the pressure force has to balance the 
gravitational force at the sonic point. 

The sonic point r 3 can therefore be obtained from the 
given potential profile and the given speed of sound by 



i>'{Ts)= 2 -^. 

I S 



(8) 



1 There are some exceptions to this line of arguments, for exam- 
ple when the force is a function of velocity. 
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Generally, the solution of this equation is expected to 
be implicit and it should therefore be solved numerically. 
Nevertheless, one can anticipate a solution of order r s ~ 
GM/c 2 (where M is some typical mass within the sonic 
radius). In other words, the escape velocity from the sonic 
radius will generally be of the order of the speed of sound. 

We next wish to obtain the velocity profile of the wind. 
By integrating the Euler equation over r, with the steady- 
state condition v = 0, we obtain 



v 2 (r) + c s lnp(r) + <j>{r) = const. 



(9) 



a constant in r. By integrating the continuity equation, with 
the steady-state condition p = 0, we obtain 



In p(r) + 2 lnr + lnv(r) = const. 

When combined, the two equations yield: 



1 v 2 (r) 
2~cT 



\nv(r) — 2 lnr + 



(r) 



= const. 



(10) 



(11) 



The constant is determined by the requirement that v = c, 
at the sonic point. In exponential form, the equation be- 
comes 

v 2 {r) 

I — | _ ; ('XI) I . — — — — 



v(r) 

— — exp i — 

c s \ 2c% J \ r J \ c% c; 

This equation relates the velocity profile and the potential 
profile, where r s is obtained from eq. (0- 

The mass loss rate can now be evaluated by equating 
its value at the sonic point and at some small intermediate 
radius ri where v <C c s , using the wind velocity profile of 
eq. (1121 . We obtain 

M = 4np s r 2 c s = A%pir 2 Vi 



(13) 



Aitpir s c a exp 



The wind solution is not strictly valid at r — > 0, but since 
we chose n such that u C c s , we can relate p(n) to the gas 
density po at r = via the static solution to eqs.0 through 

ED 



/'n = rx]> { — {■ •(;•, I - 0(0)] 
Thus, the mass loss rate is 



M = 4npor a c s exp 



- 0(0)] 



(14) 



(15) 



The total mass of gas within r s , using eq. 114H . is ap- 
proximately 

M gas {r s ) = 4np ( r 2 exp (-Ltyfj.) - 4>{0)]j dr. (16) 



The evaporation time scale, defined by r evap = 
Mgas/M, can now be obtained from eq. 1161 and eq. 115H : 

(^exp([0(r s )-0(r)]/^)dr. (17) 

Note that it does not involve the values of the potential or 
the density at the halo center. 

To obtain the above result, we implicitly assumed that 
the static solution, eq. 1141 . is valid for the density out to r s . 
This is a sensible approximation as long as the bulk of the 
mass at r < r s has v <C c s . Under this limit, where most of 



the mass is almost at rest, the total mass estimate of eq. 1161 
is valid. If we only have v ^ c s , a large fraction of the mass 
evaporates on an acoustic time scale - the sound crossing 
time of the characteristic length r s . In this case we do not 
have a steady stat e and the system shou ld be treated in 
another limit, as in lBarkana fc Loebl dl999T) . Quantitatively, 
the requirement that our solution be valid is therefore 



Tevap Tacoust 



r s /c s 



(18) 



In order to apply the above result, we need a potential 
profile and the value of c s as input. We then find the sonic 
point by eq. ©, and then use eq. I15H or eq. 1171 1. 



3 APPLICATION TO REALISTIC HALOES 
3.1 An NFW Profile 

We apply the above evaporation analysis to the typi- 
cal profile of galactic haloes, as seen in cosmological N- 
body simulations of the standard ACDM cosmology. The 
dark-matter density is assumed to follow the NFW profile 
jNavarro et al. Ill997f) : 



Pdm(r) = 



4p c 



'(1 + x) 2 



(19) 



where r c is a characteristic inner radius and p c is the density 
at r c . The mass profile is obtained by an integral over r: 



M dm (r) = lQirp c rlA{x), A(x) = ln(l + x) 



(20) 



We may replace below the parameter p c by M c , the mass 
enclosed within r c , which is given by 



87rp c rc(ln4 



(21) 



Assuming that the gravity is dominated by the dark 
matter, the gravitational force is 



fgrav{f) 



GM dm (r) 



-l6Gnr c p c ^Y- 



and the potential is: 



fgrav{r) = l&GlY p c r c B(x) 



where 



B(x) 



x — ln(l + x) 



(22) 



(23) 



(24) 



The sonic point r s is the solution of eq. (JHJ, which takes 
the form: 



A(x s ) _ (In 4- 1) _ GM c /r c 

X s 1p ' _ C 2 S 



(25) 



The dimensionless parameter ip characterizes the general be- 
havior of the gas in the given system. It relates the gravi- 
tational potential at the core radius to the specific thermal 
energy of the gas. In other words, it is the ratio squared 
between the circular velocity at the core radius and the 
isothermal speed of sound. Depending on the value of ip, 
the system could be in one of the following three regimes. 

2 One may loosely use the term "core" radius even when it is 
actually a steeper "cusp" . 
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lo gio( x ) 

Figure 1. The solution for the evaporation of gas from an NFW 
potential, with ip = (GM c /r c )/Cg = 2.5, in the regime where the 
evaporation is through a steady wind. The sonic point is at x s = 
9.18. Shown are the NFW potential <j>(x) in units of ej, the gas 
velocity in units of c s , and the gas density p(x) in units of its value 
at the sonic point. Note that since the sonic point is likely to be 
comparable to the virial radius, our spherical solution may become 
invalid outside the sonic point. This does not change the estimates 
of the mass loss, but it may change the expected gas velocity at large 
radii. 
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Figure 2. The dimensionless functions of ip for NFW haloes, where 
ip characterizes the ratio of potential to thermal energy. For ip < 1-78 
there is no steady state solution, and evaporation takes place after 
a few sound-crossing times. The function Xs(tp) is the radius of the 
sonic point in units of the core radius. The function fi(4>) is the 
gas mass within the sonic point in units of 47r/3pnrj? where po is 
the gas density at r = 0. The function /2 ("</>) is the mass-loss rate 
in the above mass units per sound crossing time r c /c s . The ratio 
/1W//2W is the thermal evaporation time scale in units of r c /c a . 



A value of ip 3> 1 implies that the gas is tightly bound and 
the evaporation is negligible. When ip > 1 but ip is not too 
large, we expect a steady-state evaporation by a thermal 
wind over many dynamical time scales. On the other hand, 
when ip ;$ 1, we expect no steady-state behavior and the 
gas evaporates rapidly over a few dynamical time scales. By 
solving eq. 11251 . we obtain the dimensionless sonic radius 
x a (ip). 

Once x s is found, we can obtain the actual wind solu- 
tion. To do so, we take the NFW potential, eq. 1231 . and plug 
it in eq. 1121 . which is an implicit equation for v(r)/c s . This 
equation can then be solved numerically. The solution for 
the wind velocity profile v(x), together with the gas density 
profile p(x) and the NFW potential, are plotted in fig.^for 
the specific case of ip — 2.5, for which x 3 — 9.18. An inter- 
esting point to note is that this sonic radius is comparable 
to the virial radius in a typical NFW halo, w hich is on the 
order of R v ~ 10r c fsee lBullock et ql~ll2001all . This implies 
that the wind solution at the sonic point may be perturbed 
by deviations from spherical symmetry due to neighboring 
haloes. This may have a noticeable effect on the velocity at- 
tained by the evaporated gas, but since the sonic point is 
determined by conditions internal to it, we expect no signif- 
icant effect on the mass-loss rate itself. This is because, by 
definition, material past the sonic point is moving too fast 
to acoustically affect conditions up stream. 



A useful dimensionless function of ip is the gas mass 
within r s in units of (4-7r/3)po?"?: 



M gas 



(47r/3)p r c 3 

.«>) 



(26) 



2 

x exp 



20 



(ln4- 



-B(x) 



dx, 



where M gas is taken from eq. 1161 with the NFW potential 
from eq. 1231 . Note that ip enters both in the integrand and 
through the upper limit of the integration x s (ip). 

Another useful dimensionless function of ip is the mass- 
loss rate in units of (47r/3)porj? per the gas acoustic time 
scale r c /c s : 



M 



— 3 exp 



(47r/3)por?/(r c /c s ) 

1 2ip 

2 ~ (In 4- 1) 



(27) 



B[xs(n 



The evaporation time scale, for an e-fold mass loss, is 
therefore, as in eq. 1171 . given by: 



7~evap — 



Mga S _ r c fijtp) 
M ~ c s f 2 (ip)' 



(28) 



The solution of x a (ip), together with the functions, 
fi(tp), h{ip), and fi{4>)/ f2(4>), are depicted in fig. [5] as a 
function of ip. The ratio fi{ip)/ f2(ip) gives the evaporation 
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time scale in units of the sound crossing time. Take, for ex- 
ample, the case where ip = 3. Most of the gas is dynamically 
bound, as the potential energy dominates over the thermal 
energy, but the evaporation process could be important if 
enough time is available. In this case, /i (VO//2 W) ^ 50, 
such that it takes about 50 sound crossing times of the core 
radius to evaporate the gas. 

3.2 Haloes of different mass and velocity 

Since we wish to compare the above results to actual cosmo- 
logical haloes, it is convenient to express the NFW profile 
parameters in terms of the "virial" quantities: the radius R v , 
the mass interior to it M v , and the circular velocity at the 
virial radius V v . Also, instead of using the more physically 
relevant acoustic time scale, we should express the evapora- 
tion time in units of the Hubble time at the corresponding 
cosmological epoch. 

At redshift z, the virial mass is related to the virial 
radius through the universal top-hat overdensity A at viri- 
alization: 

47T 



C = R v /r c 



(37) 



AI, 



-A{z)p u R v , 



(29) 



where p u is the average universal density at z. The corre- 
sponding virial velocity is defined by 

GM V 



V 

1 v 



(30) 



For the family of flat cosmologies ( f2 m + Qa = 1), the 
virial overdensity can be approximated l|Brvan fc Norman"! 
1998) by: 



A(z) 



'18tt 2 + 82[fi m (; 



39[fi m (z) - l] 5 



(31) 



Q m (jz) 

where Slm{z) is the mass density parameter at redshift z. It 
is related to today's values of the cosmological parameters 
by 

ftm0(l +Z) 3 



o m0 (i + z)3 + i-a, 



(32) 



At the high redshifts of relevance to us here, £7 m is close to 
unity, and A(z) ~ 180. For the standard ACDM cosmology, 
with Qmo = 0.3 and Qaq = 0.7, today's value is A(z = 0) ~ 
340. 

If we define a factor of order unity by 
n m0 A(z) ( h ^ J 



(&)' 



(33) 



0.3 200 

where h is the Hubble constant in units of 100 kms -1 Mpc -1 , 
we can write the above virial relations (to an accuracy better 
than 1%) as 



M s = 0.342Fi(l + Z ) 3 J R?o, 



(34) 
(35) 

M 8 = 6.06-F\" x/ '(l + z)- a/ *Vi , (36) 



V 10 = 0.383Ft /2 (l + zf /2 R 10 , 



where M 8 = M V /10 8 M Q , R w = 7? v /10kpc, and V 10 = 
K/lOkms" 1 . 

To relate the virial variables to the core characteristics 
of the NFW profile we use the concentration parameter, de- 
fined by 



At a given M v and z, we use here the average conc entra- 
tion parameter as derived bv lBullock et al. I i2001ah from 
an iV-body simulation of the ACDM cosmology. It can be 
approximated by 



C(M V , 



20(1 + 2) 



M v 



10 n M, 



(38) 



Once we combine the above relations to the results ob- 
tained for the thermal evaporation of the gas, we are in a 
position to evaluate the evaporation process of typical cos- 
mological haloes. In particular, we can compare the evapo- 
ration time scale, by which the halo has lost an e-fold of its 
gas mass, to the Hubble-expansion time scale at the onset 
of ionization, Zi on . This is depicted in fig. where we plot 
T eV ap/THubbie(zion) for a set of initial and terminal redshifts 
for the ionizing radiation. From the figure, it is apparent that 
for earlier values of Z[ on , the thermal evaporation process can 
evaporate gas from haloes with a somewhat higher virial ve- 
locity. This is mainly because at a higher redshift a halo of a 
given V v has a smaller ip, roug hly ip oc (l + z)-°- s . [It comes 
about because, for a given V v , R v oc M v oc (1 + z) _3//2 and 
roughly C oc (1 + z)" - 8 , such that r c = R v /C oc (l + z)~ - 7 . 
On the other hand, M v oc (1 + z)~ 3/2 , and M c /M v is only 
a weak function of C. Together, with c s the same at all 
redshifts, we get that roughly ip oc M c /r c oc (1 + z)~ 0,8 .] 
We see in Fig. [5] that /1 //2 is varying exponentially with 

V*> /1//2 x e b< - 1+z ^ 1 , so we expect a high reduction in 
/1//2 at a high z. The multiplication by r c oc (1 + z)~ ' 7 
in eq. (1281 makes the decrease of r e „ ap with z even steeper. 
The corresponding decrease in THubbU, which at early times 
is oc (1 + z)~ 3 / 2 , weakens the decrease of T eva p / THubble with 
z, but it is not steep enough to compensate for the exponen- 
tial decrease in Tevap- The net result is therefore an increase 
in the evaporation efficiency as a function of z. For exam- 
ple, if t = 10 4 K, the critical velocity for e-fold loss in one 
Hubble time grows from V v = 18.8 km s -1 at Zi on — 1 to 
V v — 22.0 kms^ 1 at z- lon = 20. However, if the ionizing flux 
remains effective until, say, 2 = 1, the critical velocity for 
Zion ~ 20 becomes V v = 27.7kms _1 . 

The terminal redshift, when the ionizing radiation is 
switched off, is also a factor in determining the critical halo 
virial velocity from which gas can evaporate. For example, 
if T = 10 4 K and z- lon = 10, haloes that lose an e-fold of their 
gas in either one Hubble time, by z = 5, z = 2, or Z = 1, are 
of Vv « 21.5, 22.3, 25.5, and 26.4 km s" 1 respectively. 

On the other hand, the temperature at which the gas 
can be kept heated is more important at determining the 
evaporation losses. This is because it determines the speed 
of sound. Since c s oc VT, all typical velocities in the system 
scale in the same manner. Namely, by increasing T from 
10 4 K to 2 x 10 4 K, the critical velocity for e-fold evaporation 
becomes larger by a factor of ~ v2. For example, if Zi on — 
10, the critical velocity for an e-fold loss in a Hubble time 
grows from V v = 21.5kms _1 for T = 10 4 K to V v = 29.3 for 
T — 2 x 10 4 K. 

The intersections between the acoustic curves and the 
evaporation curves plotted in fig. |S| correspond to the low- 
est velocities for which there is a consistent steady state 
thermal evaporation, and e-fold mass loss over a dynami- 
cal time. These velocities are somewhat larger than the ve- 
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V vir [km/s] 

Figure 3. The evaporation time of haloes with different virial ve- 
locities V v - The calculation is carried out at different redshifts z-, OI , 
when the gas is assumed to be photo-ionized (as marked). The evap- 
oration time is expressed in units of the Hubble time at that z- lon . The 
sets of curves on the left and on the right correspond to T = 10 4 K 
and 2 X 10 4 K respectively. The velocity at which the curve reaches 
Tevap/THubble = 1 marks the critical virial velocity of a halo that 
loses an e-fold of its gas in one Hubble time. Also marked on each 
curve are the Hubble times at later redshifts z en d at which the ion- 
ization terminates. For example, a halo which is ionized to 10 4 K at 
-Zion = 10 (the solid curve in the left set), and is subject to an effec- 
tive ionizing flux until z ond = 1 (the dot on the curve), will decrease 
its gas mass by an e-fold if its virial velocity is V v = 26.4 km s . If 
T = 2 X 10 4 K, a halo of V v = 36.8 krns" 1 will suffer a similar mass 
loss by z = 1. The lower set of curves is the "acoustic time scale" 
for the different redshifts. Our wind calculation is valid as long as 
this time scale is much smaller than r evap . 



locities corresponding t o instantaneous mass loss found by 
iBarkana fc Loebl lll999l) using binding energy arguments. 

The fraction of gas not evaporated decreases exponen- 
tially with time, such that over a given time At, this frac- 
tion is exp(— At/Tevap)- Thus, using the evaporation time 
scale, we can estimate the amount of gas that will evapo- 
rate over a Hubble time. This is described in fig. 2] which 
plots the fraction of evaporated g function of the 

virial mass of the halo, or the virial velocity in the inset. 
Also plotted is the fraction of gas which becomes unbound 
and "instantly" evaporates on a dynamic al time scale once 
ioniza tion is switched on, as calculated bv lBarkana fc Loebl 
(1999). 

From the figure, it is evident that by letting the slow 
process of thermal evaporation take place, a given relative 
mass loss is obtained in haloes which are typically 10 times 
more massive. For example, only haloes of ,6 2 x 1O 7 M0 
( lOkms -1 ) lose 70% of their gas mass from the instan- 
taneous evaporation, while the thermal wind ensures that 
haloes with < 2 x 10 8 M Q ( < 22kms _1 ) lose 70% of their 
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Figure 4. Fraction of gas mass evaporated from haloes due to 
re-ionization at z- lon = 8 as a function of halo mass or virial ve- 
locity (inset). The two pairs of curves on the right are the results 
of evaporation by thermal wind continuing until z cn( j = 4 and 1, 
for T = 10 4 and 2 X 10 4 K as marked. The two pairs of curves on 
the left are the gas fraction evaporated instantaneously according to 
IBarkana & Loebl 1 1999). In each pair, the curves are for T = 10 4 
(left) and 2 X 10 4 -ft" (right). The pair of curves on the very left, 
showing less evaporation at small masses, represents the inclusion of 
self shielding, which enables the existence of a dense core inacces- 
sible to the effects of the external ionizing radiation. We find that 
the continuous evaporation by wind increases the critical mass for 
evaporation by an order of magnitude, and it sharpens the transition 
between efficient and inefficient evaporation. 

gas mass, if T = W 4 K . For T = 2 x 10 i K, the equivalent 
mass increases by more than a factor of two. In the case 
of haloes which suffered higher mass loss fractions, we find 
that the inclusion of thermal evaporation can increase the 
equivalent cutoff virial mass by almost 100-fold. This is be- 
cause the thermal wind can even evaporate the dense cores 
which otherwise are shielded from the ionizing radiation and 
therefore do not evaporate dynamically. 

We also see in the figure that the process has a sharp 
cutoff, that is, either a halo can or it cannot retain its gas. 
Only a very narrow range of halo masses partially retain 
their original gas content. Also apparent is that the critical 
halo mass is not sensitive to the epoch when the ionization 
becomes ineffective. We learned in Fig. [3] that it is also not 
too sensitive to the epoch when the ionization starts, but it 
is quite sensitive to the temperature at the ionized state. 

3.3 Mass budget 

Next, we wish to evaluate the total fraction of gas evapo- 
rated from haloes during the ionization period. We consider 
the distribution of halo masses at Zi on and integrate over the 
fraction of gas evaporated from each of these halos, either in- 
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stantaneously or under the assumption that the ionization 
remains effective until z en d = 2. The results are shown in 
Fig. El 

The distribution of halo masses is estimated us ing the 
Sheth-Tormen modification jSheth fc Tormenl|l999T) of the 
Press-Schechter approximation. Assuming a cosmological 
model and Gaussian initial fluctuations with a given lin- 
ear power spectrum, the fraction of mass in haloes of virial 
mass larger than M is estimated by 




where v = 5 c /[D(z)a(M)}, with S c = 1.69, D(z) the lin- 
ear fluctuation growth rate, and a(M) the rms ampli- 
tude of linear density fluctuations today top-hat smoothed 
over mass M. We assume the standard flat ACDM cos- 
mology with O m p = 0.3 a nd JIao = 0.7. Following 
ICarroll. Press fc Turner! il992l) . we approximate 

_ Qtio + (3/2)n m0 On(z) u 

{) ~ n m0 [n m ( z )v + (3/2)n m (z)](i + z) ,( - 

with Q m (z) from eq. 13211 . The rms fluctuation a(M) is de- 
rived in the standard wa y using a top-hat wind ow and the 
CDM power spectrum of iBardeen et al. 1 il98rJ) with bary- 
onic density f2 D = 0.044, Hubble constant h = 0.7, power 
index n = 0.95 and the amplitude of the power spectrum 
normalized by as = 0.85 (motivated by the best fits from 
the CMB observations by the MAP satellite) . 

In Fig. El for every value of Zi on , the gas is divided 
into three different populations. Below the shaded area is 
the gas that either has evaporated from mini-haloes of 
V v < lOkms" 1 (where the inefficient cooling does not al- 
low any stars to form) or "field gas" that has not been in 
any virialized halo to begin with. The shaded area refers to 
the fraction of the gas at Z[ on that will be photo-evaporated 
by z C nd = 2 from haloes of V v > lOkms -1 . The region 
above the shaded area is the rest of the gas, bound to mas- 
sive haloes and not to be photo-evaporated. The evaporated 
fraction is further divided into gas that escapes on a dynam- 
ical time scale at Zj on versus gas that flows out continuously 
in a thermal wind until z cn d = 2. The results for the two 
possible gas temperatures are shown. 

We find that quite independently of Zi on , the amount of 
gas which collapsed into virialized objects and later evap- 
orated is about 6% if the temperature of the gas was kept 
at ~ 2 x 10 4 K, or about 4% if T was kept at ~ 10 4 K. The 
instant evaporation involved 1 to 1.5% of the gas. In other 
words, the photo-evaporation process recycled gas in dwarf 
galaxies, and it can therefore have implications on the metal- 
licity of the IGM. 

3.4 The assumption of isothermality 

The results thus far described should be considered cau- 
tiously until we verify the basic assumption behind the evap- 
oration process — that the gas is kept isothermal during the 
evaporation process. This implicitly assumes that the typi- 
cal time scale for gas evaporation is longer than the time it 
takes the expanding gas to heat and establish equilibrium 
with the radiation field. If this assumption fails, the gas will 
expand adiabatically instead and cool rapidly with distance. 




Figure 5. "Mass budget" at the onset of reionization z- lon . The 
area below the shaded area refers to gas which is either originally 
associated with haloes of V v < lOkms - 1 or not associated with 
any virialized halo. This gas is not involved in star-forming galaxies. 
The gas above the shaded area is bound to massive haloes and it will 
not photo-evaporate by z cn( j = 2. The shaded area refers to the gas 
that will evaporate from V v > lOkms" 1 haloes. The top two curves 
refer to continuous evaporation by z en d = 2 for the two possible 
temperatures as marked. The shaded area between the dashed lines 
refers to gas that escapes on a dynamical time. One may assume 
that the gas in the shaded strip is recycled in dwarf galaxies and 
returned to the IGM. 



Without the extra energy input from the radiative heating, 
the net result will be significantly less mass loss from given 
haloes. 

Under the assumption that the ionizing radiation is 
strong enough to keep the gas hot and in equilibrium be- 
tween heating and cooling, any additional assumption re- 
garding the heating rate is equivalent to an assumption re- 
garding the cooling rate. Thus, we can test instead the con- 
dition that the cooling time is short enough, as compared 
with the time it takes the gas to outflow from the potential 
well and reach the sonic point. This is typically the sound 
crossing time at the virial ra dius. 

From the results of ISuth e rland fc Dopital dl993f) : 
iKatz. Weinberg fc Hernauistl <ll996f) . we find that for tem- 
peratures in the range 10 4 K and 10 5 K and for pri- 
mordial composition the cooling rate is A = 5 x 
10" 23 A nierg s" 1 cm , where uh is the number density 
of the gas and Ao ~ 1 to within a factor of two (If the 
gas is contaminated by metals, the cooling rate could be 
even faster). The characteristic cooling time can then be 
defined by r coo i = (3/2)Csp/A, where c s is the isothermal 
speed of sound, corresponding to a temperature T. The 
sound crossing time at the virial radius defines the dy- 
namical time, Tdyn = R v /c s . We write the gas density as 
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p = Abfibp cr i ti o(l + z) 3 , where Af, is the overdensity of 
baryons at the virial radius, and p cr it,o is the critical den- 
sity of the universe today. Adding these expressions together 
results with: 

S^- 5R (d)(^) (i + ^ (41) 

where F2 is a factor of order unity, 

In the range of interest we have V v ~ 20 kms" 1 and z > 1—2, 
so Ai, > 10 would guarantee efficient cooling. We can eval- 
uate A;, i n the case where no win d is present by inspecting 
Fig. 2 of iBarkana fc Loebl l|l999T ). which shows that A& is 
typically £ 10. The presence of a wind can only increase the 
density at the sonic point, because in this case the density 
falls off at large distances as a power law, while in the no- 
wind solution it falls off exponentially. We therefore expect 
Ab > 10, and conclude that the condition for isothermality 
should be valid. Before z ~ 4, we expect the cooling time to 
be shorter than the dynamical time by more than an order 
of magnitude. 



4 DISCUSSION 

We analyzed the continuous evaporation from dwarf galaxies 
by thermal winds during the period when an ionizing flux 
was effective. We found that the critical virial velocity for 
significant evaporation is significantly higher than computed 
before assuming an instantaneous mass loss. For example, 
if the ionization starts at Zion = 10 and is maintained until 
Zend = 2, a mass loss of one e-fold occurs in haloes of V" eV a P — 
24 to 34 kms -1 , depending on whether Hydrogen or also 
Helium are reionized respectively. One e-fold of gas is lost 
within the first Hubble time at z = 10 for haloes of V C va P — 
20 to 29 kms -1 respectively. 

Before we continue with the discussion of possible impli- 
cations, we should mention the main caveats in our analysis. 
First, we assumed spherical symmetry for the system of dark 
halo and gas. This is generally not strictly valid, because 
the galaxies are typically triaxial configurations which may 
also involve significunt substructure. The fact that the equi- 
potential surfaces tend always to be rounder than the equi- 
density surfaces helps making the spherical approximation 
sensible for our analysis. Nevertheless, this is a limitation on 
the accuracy of our results. 

The other caveat, as discussed before, is the fact that 
the sonic radius obtained is typically comparable to the halo 
virial radius. On scales larger than the virial radius, one may 
expect the potential wells to deviate from the NFW profile 
and be perturbed by neighboring structure. In fact, the po- 
tential well need not be spherically symmetric there, nor 
should it even continue to increase at very large distances. 
This implies that one should not take too seriously the wind 
solution at r > r v ~ r 3 . For example, the formal solution 
shown in Fig.^indicates that the wind continues to acceler- 
ate past the sonic point to more than twice the sound speed 
at large radii. This will not necessarily be the case. Never- 
theless, once the wind turns super-sonic, it obviously cannot 
affect conditions present upstream. This implies that as long 



as the sonic point or regions internal to it are not affected by 
nearby haloes, conditions at the sonic point which determine 
the mass loss rate, will not change. 

As long as the IGM is ionized, accreting gas is not 
expected to quench the steady thermal winds. This is be- 
cause the Jeans scale corres ponds to haloes of virial velocity 
~ 30 kms" 1 or larger (e.g. IBarkana fc Loebl EoOll) . imply- 
ing that no significant gas infall is expected into haloes that 
support a steady wind. 

The upper critical velocity for evaporation is thus com- 
parable to the Jeans scale, and significantly larger than the 
cooling barrier at ~ 10 kms -1 , below which it is believed 
that luminous galaxies cannot form. This implies that every 
halo in the crude range 10 — 30 kms -1 is likely to be af- 
fected by the photo-evaporation process. The mass fraction 
in such haloes is on the order of 5% at any given time after 
z ~ 10. Some of them eventually merge into bigger haloes, 
but some end up as dwarf haloes today. Among these un- 
merged haloes, those that managed to form some stars prior 
to Zion possibly make low-surface-brightness dwarf galax- 
ies which are dominated by an old population of stars and 
with no original gas retained. Young stars can form from 
newly accreted gas after the ionization stops. These objects 
may resemble the dwarf spheroidals of the Local Group. 
The other haloes in the range 10 — 30 kms -1 , those that 
are dynamically younger and did not form stars prior to 
2:i on , were doomed to never form stars, at least until z C nd- 
They are likely to end up as completely dark haloes. These 
are therefore not limited to the mini haloes below the cool- 
ing barrier of ~ 10 kms -1 , but they rather extend up to 
~ 20 — 30 kms -1 . This may help explain the apparent 
discrepancy between the predicted mass function and the 
observed luminosity function at the faint end. While cer- 
tain variants of the currently accepted cosmological model 
may predict a number density for subhalo satellites in big 
haloes that matches the observed population of Local-Group 
dwarf galaxies in the range V > 20 - 30 kms -1 , it is much 
harder for such models to recover the small number at 
V < 20kms -1 dZetner fc Bullock! |2003> . This can be ex- 
plained by the very efficient photo-evaporation of gas tak- 
ing place in haloes below I4 V a P ~ 20 kms -1 , as found by our 
wind analysis. 

Given our results, we may attempt to address addi- 
tional aspects of the possible scenario for dwarf spheroidal 
galaxies. They are constrained to form in haloes which lie 
in the crude range 10 < V < 30 kms -1 . The gas in haloes 
of V < 10 kms -1 cannot cool to form stars at any early 
epoch because the cooling rate has a sharp drop below 
T ~ 10 4 K; any molecular hydrogen that is necessary for 
cooling at T < 10 4 K is likely to be dissociated even be- 
fore reionization. Efficient cooling by Hydrogen recombina- 
tion at T £ 10 4 K can lead to an early burst of stars in 
haloes of V Si 10 kms -1 . The associated supernovae blow 
out part of the gas, as they do with a grad ually decreas- 
ing efficiency in all haloe s of V < 100 kms -1 (iDekel fc Silkl 
ll986tlDekel fc Wool2003l) . In V < 30 kms -1 haloes, the rest 
of the gas photo-evaporates in a few Hubble times, as com- 
puted above. No further gas can fall in as long as the ion- 
izing flux is effective, which may last until z ~ 1 — 2. This 
leads to a gas-poor, old-population stellar system. Any re- 
cent star formation in these dwarf spheroidals must be due 
to gas falling in after the ionizing flux became ineffective. 
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An inspection of the global properties of dwarf galax- 
ies indicates that the dwarfs in the range 10 — SOkms" 1 
span more than two orders of magn itude in luminosity, or 
stellar mass M» foekel fe Wodl2003l) . They show no signif- 
icant correlation between velocity and luminosity, contrary 
to the irregular, disky dwarfs that are typically larger than 
~ 30kms _1 and tend to lie on an extension of the tight 
Tully-Fisher relation defined by the brighter galaxies. This 
observed property is consistent with the scenario sketched 
above for the dwarf spheroidals. The spread in A/* within 
this family represents variations in quantities other than 
halo velocity, such as the proximity of the earlier star for- 
mation epoch to Zi n- The spread in M* for a given halo V 
becomes much smaller for the larger dwarf irregulars and 
bright galaxies because haloes above V eva , p retain a larger- 
fraction of their gas and being above the Jeans scale they 
allow further accretion of more gas. Star formation is thus 
allowed to occur in these objects also after z lon , such that a 
large fraction of the remaining gas can turn into stars. This 
brings the value of M* close to the value predicted by the 
energy requirement for gas removal by supernova, which is 
a strong function of the halo virial velocity: M»/M oc V 2 
JPekel fc Wooll2003l) . 

The above scenario may also explain why the dwarfs at 
the low end, 10 — 30kms _1 , are typically spheroidals with 
low spin, while the larger dwarfs tend to be disks. Accord- 
ing to the above scenario, the stars in the smaller galaxies 
are made from the gas that cooled first, namely the gas 
residing in the inner halo. This gas is typically of specific 
angular m omentum significantly lower than the average for 
that halo {Bullock et al. 1 12001 ri iMaller fc Dekel l2002h . In 
fact, the scenario predicts a correlation between M* and 
stellar spin parameter because the low-M* dwarfs typically 
form later, closer to Zj on . The smaller fraction of gas involved 
in the fainter dwarfs must have originated from inner halo 
regions which host the low-end tail of the specific angular- 
momentum distribution. There is a hint for such a corre- 
lation in the Local Group dwarf spheroidals JPekel fc Wool 
120031) . 

The photo-evaporation process discussed in this paper 
may involve a non-negligible fraction of the gas in the uni- 
verse. If it evaporated from dwarfs after they formed stars, 
it could carry metals with it and enrich the IGM. This could 
provide a lower limit to the actual enrichment, which is prob- 
ably dominated by supernova blowout. We estimated for the 
standard ACDM cosmology that at any given time about 
5% of the mass is in haloes that can allow significant evap- 
oration. Integrated over time, we find that by z ~ 1 about 
40% of the mass is in haloes of V > lOkms - . A signifi- 
cant fraction of this mass must have been in haloes in the 
evaporating range for more than a dynamical time sometime 
during the merger history. This implies that up to 40% of 
the gas has been evaporated from dwarf galaxies, and could 
carry metals. 
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